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ABSTRACT 



f 



Aims. We construct a theoretical model to predict the number of orphan afterglows (OA) from gamma-ray bursts (GRBs) 

triggered by primordial metal-free (Pop III) stars expected to be observed by the Gaia mission. In particular, we consider 

primordial metal-free stars that were affected by radiation from other stars (Pop III. 2) as a possible target. 

Methods. We use a semi-analytical approach that includes all relevant feedback effects to construct cosmic star formation 

history and its connection with the cumulative number of GRBs. The OA events are generated using the Monte Carlo method, 

and realistic simulations of Gala's scanning law are performed to derive the observation probability expectation. 

Results. We show that Gaia can observe up to 2.28 ± 0.88 off-axis afterglows and 2.78 ± 1.41 on-axis during the five-year 

nominal mission. This implies that a nonnegligible percentage of afterglows that may be observed by Gaia (~ 10%) could have 

Pop III stars as progenitors. 
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1. Introduction 

The first stars (hereafter. Pop Ill-primordial metal-free) 
in the Universe are thought to have played a crucial role in 
the early cosmic evolution by emitti ng the first light an d 



producing the first heavy elements (jBromm et al. I I2OO9I) . 



Understanding such objects is very important since their 
detection would permit the pristine regions of the Universe 
to be probed. However, there has been no direct observa- 
tion of the so-called Pop III stars up to now. 

Pop III stars may produce coUapsar gamma-ray 
bursts (GRBs) whose total isotropic energy c ould be 
Rri 2 orders of magnitude larger than average ( Barkov 
2010l:lKomissarov fc Barkovll2oiol:lMeszaros fc Reedl20ld 



Suwa fc Iokall201ll : lToma et al.ll201l[) . Even if the Pod III 

star has a supergiant hydrogen envelope, the GRB jet can 
break out of it because of the long-lasting accretion of th e 
envelope itself (jNagakura et al.ll2012l : ISuwa fc Iokall2011l ). 
It is of great importance to study the rate and detectabil- 
ity of Pop III GRB prompt emissions and afterglows in 
current and future surveys. We explore here the possi- 
bility to observe these objects through their afterglows 



(jToma et al.ll201l[ ). Observations of GRB afterglows make 
it possible to derive physical properties of the explosion 
mechanism and the circumburst medium. It is intriguing 
to search for signatures of metal-poor stars in the GRB 
afterglows at low and high redshifts. 

GRB optical afterglows are one of the possible tran- 
sients to be detected by the Gai£0 mission. Recently 



Japeli fc Gombod (120111 ) explored the detcctability of 



such afterglows with Gaia using a Monte Carlo approach 
that inspired us. As the GRB jet sweeps the interstellar 
medium, the Lorentz factor of the jet is decelerated and 
the jet starts to expand sideways, eventually becoming 
detectable by off-axis observers. These afterglows are not 
associated with the prom pt GRB emission and are called 
orpha n afterglows (OA) (jNakar et al. 2002 : Rossi et al 



20081) . 



De Souza et al. (2011) showed that, considering 
existH specifications, we can expect to observe a max- 
imum of « 0.08 GRBs with z > 10 per year originat- 
ing from primordial metal-free stars (Pop III.l) and « 20 
GRBs with 2 > 6 per year coming from primordial mctal- 
frec stars that were affected by the radiation from other 
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stars (Pop III. 2). In the context of the current Swif^satel- 
hte, ~ 0.2 GRBs with z > 6 per year from Pop III. 2 stars 
are expected. These numbers reflect the fact that, com- 
pared to Pop III.l stars, Pop III. 2 stars are more abun- 
dant and can be observed in a lower redshift range, which 
makes them more suitable targets. In the light of such 
results, the calculations presented here will focus on Pop 
III. 2 stars alone. 

Searches have been made of OAs by both X-ray surveys 
Grindlavl 1999t Greiner et al. 2000 ) and optical searches 



Becker et al.l |2004 iRvkoff et al.l I200I: iRau et all [200I 



Malacrino et allbOOTlT The purpose of the present paper 



is to calculate the Pop III. 2 GRB OA rate that might be 
detected by the Gaia mission (for more details ab out Gaia, 
see, e.g.. iPerrvman et al.ll200ll : iLindegren I2OO9I) . 

The Gaia mission is one of the most ambitious projects 
of modern astronomy. It aims to create a very precise tridi- 
mensional, dynamical, and chemical census of our Galaxy 
from astrometric, spectrophotometric, and spectroscopic 
data. In order to do this, the Gaia satellite will perform 
observations of the entire sky in a continuous scanning cre- 
ated from the coupling of rotations and precession move- 
ments called the scanning law. For point sources, these 
observations will be unbiased and the data of all the ob- 
jects bellow a certain limiting magnitude (G— 20) will be 
transferred to the ground. Certainly, galactic and extra- 
galactic sources will be among those objects. 

Typic ally, Pop III. 2 stars are formed in an initia.lly ion - 
ized gas ( Johnson fc Brornml 2006 : Yoshida et al. 2007 ). 
They are thought to be less massive than Pop III.l stars 
but still massive enough to produce GRBs. Recent re 



suits from iGreif et al.l (|201l[ ) show that, instead of form- 
ing a single object, the gas in mini-halos fragments vigor- 
ously into a number of protostars with a range of different 
masses. It is not clear up to now how this initial range of 
mass will be mapped into the final mass function of Pop 
III stars. The most likely conclusion is that Pop III stars 
are less likely to reach masses in excess of ^ 14OM0, which 
consequentl y affect the es t imate d number of GRBs from 
Pop III.l. iHosokawa et al.l ( 2011 ). performing state-of-the- 
art radiation-hydrodynamics simulations, showed that the 
typical mass of Pop III stars could be ^ A3Mq. Here we as- 
sume that this will not affect significantly the mass range 
assumed for Pop III.2 (- 10 - lOOMg). 

The paper is organized as follows. In Sect. 2, we cal- 
culate the formation rate of primordial GRBs. In Sect. 
3, wc calculate the OA light curves and their redshift 
distribution. In Sect. 4, wc discuss the details of the 
Gaia mission and derive the probability of a given event 
to be observed by Gaia. In Sect. 5, we discuss the re- 
sults, and finally, in Sect. 6, we give our concluding re- 
marks. Throughout the paper, we adopt the standard 
A cold dark matter model with t he be st-fit cosmologi- 
cal parameters from Jarosik et al. ( 201lh (WMAP-YrTp), 

= 0.267, nA = 0.734, and Ho = 71km s-^Mpc'^ 



^ http:/ /swift. gsfc.nasa.gov/docs/swift/swiftsc. html 
* http:/ /lambda.gsfc. nasa.gov/product/map/current/ 



2. GRB redshift distribution 

To estimate the formation rate of GRB s from Pop III 



stars a t a given redshift, we closely follow Ide Souza et al 
(|2011l ). Since long GRBs are expected to follow the 



death of very massive stars, their rate could provide a 
usefu l probe for c osmic star forination history (SFH) 
e.g. lTotanilll997t ICiardi fc Loebll2000l: Ifiromm fc Loebl 



2002! IConselice et al.l 120051: ICampisi et all I2OIOI l2011al: 



Ishida et al.ll2011l:lde Souza et al.ll2011l: [Robertson fc EUid 
2012i ). However, the connection between the star forma- 
tion rate (SFR) density and GRB rate is not clearly under 



stood and can be redshift dependent (e.g., lYiiksel et al 



20081: iKistler et aLll2009t iRobertson fc Ellisll2012[) . Since 
host galaxies of long-duration GRBs are often observed 
to be metal poor. Several studies connect the ori- 
gin of long G RBs with the metallicity of thei r pro 



genitors (e.g., Meszarod 
Salvaterra fc Chincarinil 



2006 



Wooslev fc Blooipl 



2007: Salvaterra et al.l 



2006 



2009; 



Campisi et al.l 1201 Ibl ). Consequently, the GRB- SFR con- 



nection could be dependent on the cosmic metallicity evo- 
lution. However, this connection is not yet completely un- 
derstood, since there is also evidence of regions within 
GRB host galaxies kn own to possess higher metallicities 
( Levesque et al. 2010() . 

Despite such uncertainties, we expect the connection 
between SFR and GRBs to be less affected by this effect 
because Pop HI stars and their environment are metal 
poor. In other words, Pop HI stars are more likely to pro- 
duce GRBs than ordinary stars. It is important to keep in 
mind that any prediction will be convolved with system- 
atic effects that w e are not tak i ng int o account. However, 
as pointed out in llshida et al.l (j201lh . the assumption is 
good enough to agree with available observational data. 

We implicitly assume that the formation rate of long 
GRBs (duration longer than 2 sec) follows closely the SFH. 
The number of GRBs per comoving volume per time can 
be expressed as 



5- 



GRB 



(z) = ?7GRB**(2:), 



(1) 



where jycRB is the GRB formation efficiency and is 
the SFR. Over a particular time interval, A^obs, in the 
observer rest frame, the number of GRBs originating be- 
tween redshifts z and z + dz is 



dN, 



GRB 



GRB 



{zy 



1 + z dz' 



(2) 



where dV/ dz is the comoving volume element per redshift 
unit. 

2.1. Star Formation History 

To estimate the SFR at early epochs, we assume that 
stars are formed in c ollapsed dark matter halos (for 



more details, please see Ide Souza et al.ll201ll ). The num 



ber o f collapsed objects is given by the halo mass func 



tion (iHernguist fc Springell l2003t iGreif fc Bromml 12006 



Trenti fc Stiavellil l2009l) . In what follows, we adopt the 
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Sheth-Tormcn function, fsr (|Sheth fc Tormenlll999l ). To 
estimate the fraction of mass inside each halo capable of 
collapsing and forming stars, we include the following im- 
portant feedback mechanisms: 

1. H2 Photodissociation 

Hydrogen molecules (H2) are the primary coolant in 
the gas within small-mass "mini-halos." H2 are also 
sensitive to ultra-violet radiation in the Lyman- Werner 
(LW) bands and can easily be suppressed by it. We 
model the dissociation effect by setting the minimum 
mass for halos that are able to host Pop III stars 
( Yoshida et al.ll2003l ). 

2. Rcionization 

Inside growing Hii regions, the gas is highly ionized 
and the temperature is ~ 10'* K. The volume-filling 
factor of ionized regions, Qmiiz), determines when the 
formation of Pop III.l stars is terminated and switches 
to Pop IIL2. To calcu late Qhtt(-z), we closely fol low 
Wvithe fc Loebl (|2003l ) as in Ide Souza et all (|201ll) . 



Metal Enrichment 

Metal enrichment in the intergalactic medium (IGM) 
determines when the formation of primordial stars is 
terminated (locally) and switches from the Pop III 
mode to a more conventional mode of star formation. 
We assume that star-forming halos launch a wind of 
metal-enriched gas at 2; > 20. Then we follow the 
metal-enriched wind propagation outward from a cen- 
tral galaxy with a given velocity Vwmd, traveling over 
a comoving distance Rwind- We estimate the ratio 
of gas mass enriched by the wind to the total gas 
mass in each halo, and then we evaluate the aver- 
age metallicity over cosmic scales as a function of red- 
shift. We effectively assu me that the so-called criti- 
cal metallicity is very low (ISchneider et al.ll2002ll2003t 



Bromm fc Loebll2003HOmukai et al.ll2005l : lFrebel et al 



2OO7I: iBelczvnski et al.ll2010[) . Therefore, Pop III stars 
are not formed in a metal-enriched region, regardless 
of the actual metallicity. 



Rollinde et al.l ( 20091 ) investigated the role of Pop 
III stars in the cosmic metallicity evolution, in particu- 
lar, the local metallicity function of the Galactic halo. 
They show that Pop III SFR should not be larger than 
3 X MQyr~^ Mpc~'^ at any redshift. We also include 
this additional constraint as an upper limit for our opti- 
mistic model. 

The top panel of Fig. [T] shows the up p er lim it for 
Pop III.2 SFR, based on Tde Souza et alj (|201lh with 
the additional constraints cited above. The Pop III. 2 
SFR is compa red with a compilation o f independent 
measures from Hopkins fc Beacom ( 20061 ) up to z ~ 6 
and from observations of color- s elected Lyman Break 



Galaxies ( Mannucci et al.l 12007: Bouwens et al. 12008 , 



20111) . UV-hIR me asurements (Reddv et al.M2008l) and 



GRB observations (jCharv et al 



2007t lYiiksel et al.ll2008t 



Wang fc Daj l2009l) at higher z ( in the figure, these will 



be refereed to as H2006, M2007, B2008, B2011, R2008, 
C2007, Y2008, and W2009, respectively). 
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Fig. 1. Top: optimistic model for Pop HI. 2 star formation 
rate (SFR) assuming a high star formation efficiency and 
low chemical enrichment. The light points are independent 
SFR determinations compiled from the literature. 
Bottom: intrinsic GRB rate diVGRB/dz, i.e., the number of 
GRBs per year on the sky (on-axis + off-axis) according to 
Eq. ([2]). This represents our optimistic model assuming a 
high star formation efficiency for Pop HI. 2, slow chemical 
enrichment, GRB formation efficiency of Jgrb = 0.001 
and a Salpctcr IMF. 



2.2. Initial Mass Function and GRB Formation 
Efficiency 

The stellar initial mass function (IMF) is critically impor- 
tant to determine the Pop HI GRB rate. The IMF deter- 
mines the fraction of stars with minimum mass that ar e 

25Mo (jBromm fc Loebl I2OO6I ). 



able to trigger GRBs, 
The /grb factor gives the fraction of stars in this range 
of mass that will produce GRBs. 

The GRB formation efficiency factor per stellar mass 

is 

JMgr 



Vgrb — /grb - 



(j){m)dm 



^^j"^ m(j){m)dm 



(3) 



where 0(m) is the stellar IMF for which we considered 
a power law with the standard Salpeter slope 0(m) cx 
^-2.35^ Miow and Mup are the minimum and maximum 
mass for a given stellar type (respectively IOMq and 
~ IOOMq for Pop HI. 2), and A'/grb is the minimum 
mass able to trigger GRBs, which we set to be 25Mq 
(jBromm fc Loebll2006l ). 

De Souza et al. (2011) placed upper limits on the in- 
trinsic GRB rate (including the off-axis GRB). In the fol- 
lowing, we set Jgrb = 0.001 and ?7grb//grb ~ 1/87Mq ^ 
as an optimistic case, consistent with their results. The 
bottom panel of Fig. [T] shows the optimistic case for in- 
trinsic GRB rate adopted in this work. 
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3. Number of Observed Orphans 

3.1. Afterglow Model 

To calculate the afterglow light curves of Pop III GRBs , 
we fo llow t he standar d presc ription from lSari et alJ (| 19981 . 
19991) and Meszaro^ ( 20061 ) . The spectrum consists of 
power-law segments linked by critical break frequencies. 
These are ly^ (the self-absorption frequency), ly^n (the peak 
of injection frequency), and i^c (the cooling frequency), 
given by 



oc 



J/c OC 



(1 
(1 
(1 
(1 



Z) 



1/2,-3/2, 



-l.-l.l/5„3/5^1/5 



z)t 



l/2„l/2 



(4) 



where g(j)) = (p — 2)/(p — 1) is a function of the energy 
spectrum index of electrons {N{-fe)d'^e oc j~Pd-fe, where 
7e is the electro n Lorentz facto r) , eg and £b are the ef- 
ficiency factors ( Meszaro^l2006l ). Eiso is the isotropic ki- 
netic energy, n is the density of the medium, and F^^tnuK 
is the observed peak flux at luminosity distance from 
the source. 

There are two types of spectra. If i/y^ < i^c, we call it 
the slow cooling case. The flux at the observer, F,^, is given 

by 



F, 



{ve^/v^Y^^ivlv^fF^^^^^, 
ivlv Yl'^F 

ivlv l-(P-l)/2p 

/ i^in J -t I/, max; 

(i^c/i^„.)"^^-''/'(^^/^^c)-^/'F,,, 



V > Vc. 

(5) 

For > j^c, called the fast cooling case, the spectrum 



IS 



F,, 



{v.ly.fl\vlv.YF,^ 

{v/v,)^'^F,, 
(ylv^r^l^F,^ 

{vJv.r^l^iyly.^)-^l''F,. 



',max5 
',max7 



max 1 ^ ^ ^m • 



(6) 

Initially the jet propagates as if it were spherical with 
an equivalent isotropic energy of i?truo = ^j^-iso/S, where 
9j is the half-opening angle of the jet. Even if the prompt 
emission is highly coUimated, the Lorentz factor drops 
7d < ^7"'^ around the time 



2.14 



5 X 1054 



\ 1/3 / ^. \ 8/3 

j (ojj ""'^'(1+^) days, (7) 



and the jet starts to expand sideways ( loka fc Meszaro^ 
20051) . Consequently, the jet becomes detectable by the off- 
axis observers. These afterglows are not associated with 
the prompt GRB emission. 

Due to relativistic beaming, an observer located at 
9obs, outside the initial opening angle of the jet {0obs > (^j), 
will observe the afterglow emission only at i i^, when 
7d 



^7' 



100 



10 



i 0.1 



0.01 



0.001 



10- 







z= 3 


Sobs = 0.0 






e„bs = 0.05 
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Fig. 2. Example of afterglow light curve at z = 3 as a 
function of observed angle, Oobs- Wc show the evolution 
of afterglow flux F{mJy) as a function of time t (days) 
and observed angle Oobs for typical parameters: isotropic 
kinetic energy Eiso = lO^'' erg, electron spectral index 
p = 2.5, plasma parameters Eg = 0.1, cb = 0.01, half- 
opening angle jet 9j = 0.1, interstellar medium density 
n = \cm~^ and frequency i/ = 5 x lO^^Hz. The horizontal 
dotted line is the integrated Gaia flux limit; solid black 
line, Oobs = 0; dashed blue line, Oobs = 0.05; dot-dashed 
red hne, Oobs = 0.1; dotted green line, Oobs = 0.20. 



The received afterglow flux by an ofl-axis observer in 
the point source approximation, valid for Op b^ ^ Oj, is re - 
lated to that seen by an on-axis observer bv ([Granot et a" 



2002 : Totani fc PanaitescrJ I2OO2I: I.Tapeli fc Gombodl2011[) 



where 



F,{Oobs,t)^eFu/d^,£.t), 



e=(l-/3)/(l-/3cos0o,,). 



(8) 



(9) 



and P ~ \J^~ V7d- '^be time evolution of the Lorentz 
factor in given by 



where 



7d(i) = 



the 



-3/8 
-1/2 



(10) 



jet br eak time 
days (|Sari et al. 



« 0.7(1 + 
1999t ). Figure 



z)(£;5i/n)i/3(0^./O.l)- 
[5] shows four examples of afterglows as a function of 
observed angle Oobs for the case of 0j = 0.1 at z = 3 for 
typical parameters described in the flgure. The flux is 
calculated for an observational frequency = 5 x lO^^Hz 
within the Gaia bandwidth. Depending on the parameters 
of the afterglow, the light curve can appear above the 
Gaia observational limits. Due to the large quantity of 
free parameters, a Monte Carlo approach is essential to 
explore the detectability of a large number of events and 
will be explained in the next section. 
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SMC 
LMC 
MW 



(Mm-') 

Fig. 3. Extinction laws usually adopted in literature: the 
Small Magellanie Cloud (SMC) law (blue line), the Large 
Magellanic Cloud (LMC) law (red dashed line) and the 
Milky Way (MW) law (green dot-dashed line). 



3.2. Dust Extinction 

A fraction of GRBs with X-ray or radio afterglows can be 
hidden by dust absorption from their host galaxies. The 
observed flux after extinction correction can be simply 
written as (see, e.g, Eliasdottir et al. 20091 ) 



f: 



dust 



-0.4A. 



(11) 



where A\ is the extragalactic extinction along the line of 
sight, A\;ext, as a function of the wavelength A plus the 
extinction from the Milk Way, A\-^mw- 

For A\-^ext, we adopted a simple Small Magellanic 
Cloud (SMC) type extinction model. The SMC model was 
already shown to provide go od fits for several GRB after- 
glows observations (see, e.g, Eliasdottir et al. I l2009h . For 
A \:MW, we use t he av erage value 0.15 from observations 
of Schadv et al.l ( 2012 ) and adopt a typical value of 0.3 



for Ay . In Fig. |21 we show the SMC extinction curve in 
comparison with other popular models. Large Magellanic 
Cloud (LMC) and Milky Way (MW). The model choice 
has no significant effect on our results, since all of them 
have a similar trend in the G band range. In Fig. 21 we 
show the effect of dust extinction in the spectral energy 
distribution (SED) of GRBs. The effect is significant in the 
G band (~ 5 x 10^"' Hz), which will considerably decrease 
the detection rate of optical GRBs, mainly at high-z. 

3.3. IGM and DLA absorption 

For high-z GRBs, much of the optical and near- 
infrared light will be absorbed by the Lya forest, 
whic h provides a powerful tool to probe the reionization 
era (iMiralda-Escudd HqqI iMesineer fc Furlanettol l2008t 



Ciardi et al.]l201l[) . The GRB 050904 at z = 6.3 was the 
first attempt to probe the IGM through GRBs at the 
epoch of reionization by using the da mping wing at wave- 
lengths redward of the Lyman break ( Totani et al. 20061 




1x10" 2x10'-' 



5x10'-' 1x10'* 2x10'* 
v(Hz) 



5x10'* lxlO'5 



Fig. 4. Example of a spectral energy distribution for ob- 
served GRBs with (dashed lines) and without (solid lines) 
extinction, assuming the SMC extinction law. We show 
the afterglow flux F{mJy) as a function of frequency 
V (Hz) for typical parameters: isotropic kinetic energy 
Eiso = lO^"' erg, electron spectral index p — 2.5, plasma 
parameters Ce = 0.1, es = 0.01, half-opening angle jet 
6j ~ 0.1, and interstellar medium density n ~ lcm~^. 
The black line represents 9obs = 0, t = 0.5 days, blue line, 
dobs = 0, t = 5 days, red line, 9obs = 0.05, t = 0.5 days, 
and orange line, Oobs = 0.05, t = 5 days. 



iKawaiet al.l[2006h . The absorption by the ne utral IGM 
can be approximated by (jMcQuinn et al.ll2008l ) 



TiGM ~ 900 kms xhi 

H{zhost)Rb 
(1 + Zhost) 



1 + Zhost 

8 

V- — V, 



3/2 



(12) 



where Va is the rest frame of the Lya line, Rb repre- 
sents the size of an HII region surrounded by an IGM 
with neutral fraction Xfji, Zhost is the redshift of the 
GRB host galaxy, and H(z) is the Hubble parameter for 
a ACDM cosmology. To estimate Xm, we use the pre- 
scription detailed in Ide Souza et al. (1201 ll) (see Fig.l). 
The optical depth of the damped Lya absorber (DLA), 
Tola = ^hio-qK;,s(1 + ^host 

)] , can be computed by 



TOLA = 7.26 



1 + Zhost 



1 



Zobs 



f AX 



\20k 



X 

-2 



(13) 



e.g., iBarkana fc Loebl 20041 ) where Vobs = c/Xobs, (1 + 



Zobs) = Xobs/Xa, Nm is the total column density of HI in 
the host galaxy, and Xobs = AA + Xa{l + z). iVni is ran- 
domly chosen assuming a cumulative distribution func- 
tion scaling as N^ \■^^, between 10^^ 



and lO^i-^c 



' (see 

IChen et alll2007l: iMcQuinn et al.ll2008l) . For each event, 
Rb is chosen from a lognormal distribution between 1- 
100 Mpc motivated by a visual inspection in Fig. 5 from 
McQuinn etldl (|2008() . 
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Fig. 5. Redshift PDF. Probability of a given event ap- 
pearing in a certain range of redshift. 



Fig. 6. Half-opening angle jet PDF. Probability of a given 
GRB to have a particular 9j. 



3.4. Mock sample 

The mock sample is generated by the Monte Carlo 
method assuming different probability distribution func- 
tions (PDF) for each quantity as explained below. The 
medium density n is randomly chosen from a flat distri- 
bution within 0.1 — lcm~^. 



3.4.1. Redshift PDF 

We generate the GRB events randomly in redshift with a 
PDF given by Eq. ©. The probability of a given GRB 
appearing at redshift z is 



J^{dNGB.B/dz)dz 



(14) 



The PDF was generated by 10^ random realizations based 
on Eqs. ^ and (fH)) . Figure [5] shows the probability of 
finding a GRB at a given redshift, indicating that a 50% 
probability of having a GRB from a Pop III star is ob- 
tained in the redshift range z ^ 7 — 11 and a 95% prob- 
ability in the range z ~ 4 — 15. Due to absorption from 
dust, IGM and DLA, the GRBs available for observation 
by Gaia are restricted to the range z ^ 3 — 7. This re- 
sults in approximately 10^ GRBs during the entire Gaia 
nominal mission, which is the value adopted in this work. 

3.4.2. Half-opening angle PDF 



Using 



an 



Yonetoku et al 



em pirica l 
( 2005h 



opening angle estimator, 
derived the opening angle 
PDF of GRBs. Their PDF can be fitted by a power-lay 
9~^. Their results seem also compatible with the universal 



structured jet model (jPerna et al.ll2003l ) . The jet opening 



angle usually rang es between 1° - 10° (iFrail et al.ll2001l 



Cenko et al. 200' 



power law in the range 
determine the PDF of 9 



For simplicity, we assume a similar 
0.01 and 



171171 
3 



imax ^ Q_2 to 



Pf 



9,(8) 



cx 



(15) 



Figure [5] shows the PDF of Oj generated by 10^ realiza- 
tions based on Eq. (jl5p . The realizations were performed 
within the range 9j = 0.01 — 0.2. The observational an- 
gle 9obs was randomly chosen between — tt. From the 
relation iJtruo = 6jEiso/2, we assume two fixed values for 
Etrue = (2.5 — 5) X 10^" ergs, which imposes the limits 
Eiso < 5 X lO^'' — lO^^ergs respectively. 



4. The Gaia mission 

The Gaia satellite will perform observations of the entire 
sky, using a continuous scanning formed by the coupling 
of rotation and precession movements, the scanning law. 
This law guarantees that each point in the sky will be 
observed several times during the mission, as it can be 
seen in Fig. [71 

Similar to what happens with CCD meridian circles, in 
the referential of the satellite's focal plane, the sky contin- 
uously moves from one side to the other while the satel- 
lite spins. During this time, the CCD charges are syn- 
chronously transferred to compensate for the sky's appar- 
ent motion and allow the integration. 

This continuous observation strategy requires an 
equally continuous reading of the CCDs. Also, since Gaia's 
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Fig. 7. Number of times each region of the sky (in galactic 
coordinates) will be observed by the Gaia satellite during 
the entire mission. 

focal plane comprises 106 individual detectorf[f|, it is not 
possible to transfer the entire content of the focal plane 
to the Earth due to bandwidth limits. So, a continuous 
analysis of the focal plane observations is also performed 
on-board, aimed at the detection of astronomical sources. 
When a source is detected, a rectangular "window" com- 
prising a few arcseconds around the detected source is cre- 
ated (its exact size and pixel binning depend on the focal 
plane's CCD column). These windows are then transferred 
to the Earth. 

For point sources, these observations will be unbiased 
and the data from all objects in the sky, bellow a certain 
limiting magnitude, will be sent to the ground. Certainly, 
among all those objects, not only galactic sources will be 
present, but also extragalactic ones. In particular, it is 
expected that point sources up to magnitude 20, in the 
Gaia passband gH, will be "windowed" and transferre43. 

As seen in Fig. [5J some of the OA events are expected 
to remain above this limiting magnitude for a certain 
amount of time. The question that remains is whether 
their duration (at G^20) is enough for them to be ob- 
served at a reasonable rate. Only two quantities play an 
important role in estimating the probability of Gaia ob- 
serving single event from a Pop III. 2: the time whichen the 
OA remains brighter than G=20, At, and the coordinates 
(Igahbgai) whcrc thc event takes place in the sky. Since 
those quantities are continuous distributions, it is neces- 
sary to analyze how the observation probability depends 



^ Fo r a diagram of Gala's focal plane, see, e.g., Ijordl et al.l 
l|2010D . 

® This Is a broad passband, which covers from 3 30-1000 nm. 
The n ominal transmission curve can be found at IJordl et al.l 
l|2010l ^. 

^ After the mission (and during the mission for some prob- 
lematic cases) , It win be possible to reconstruct a deeper image 
around each detected source. In those reconstructed Images, It 
win be possible to reach deeper magnitudes, albeit with some 
contamination from reconstruction artifacts. 



on them by building P{At, Igai, bgai)- In the present work, 
we proceed as follows. 

For a given coordinate in the sky, we start by com- 
puting the inverse Gaia scanning law to derive a transit 
time list comprising the instants when Gala's telescopes 
will be pointing at that coordinate. To be as realistic as 
possible, we adopt the Gaia Data Processing and Analysis 
Consortium's nominal implementation of it. Then, we ran- 
domly select a point in time during the entire mission life- 
time in order to place an event of a certain duration A^. 
Using the transit time list, we check if that event was ob- 
served, considering a time window of 4.4 seconds around 
each transit; this is the time needed for the signal to cross 
the detection CCD and enter the confirmation CCD. If 
there is a superposition between the event duration and 
this time window, the event is considered detected. This 
procedure is then repeated until the estimation of the 
detection probability, which is derived by simply divid- 
ing the number of detected events by the total, does not 
vary more than 1% between iterations. Finally, the whole 
procedure is repeated for each event duration At. As a 
consequence, we obtain an adequate time-sampling of the 
P{At,lgaubgai) distribution. 

For the determination of the number of OA events ob- 
served by Gaia on the entire sky, the coordinate depen- 
dency can be averaged out, allowing P{At,lgai,bgai) ^ 
P{At) ± e. This is possible because the scanning law is 
mostly known, so we can reasonably assume that the OA 
events take place randomly in the sphere. 

The procedure described above was repeated for sev- 
eral positions on the sphere, and the mean and the stan- 
dard deviation at each event duration were computed. 
To allow a good spatial sampling for the estimation of 
P{At) ± e, we tessellate the cel estial sphere at the hierar- 
chical triangular mesh, level 4 (jKunszt et al.ll2001l) . This 
means that the simulations were performed at the center 
of 2048 triangles of approximately equal areas. 

Finally, to obtain the probabilities for the whole sky, 
an additional effect must be taken into account: the struc- 
ture of our own Galaxy. Since the OAs are extragalactic 
events, the probability of observation at the galactic plane 
or bulge should be null or very small, due to extinction 
and crowding. In this work, we conservatively assumed 
a null value for the probability of OAs being observed 
at such regions of the sky (defined here as |6| < 15° for 
345° < ; < 15° and |6| < 5° otherwise). 

The final results, representing the behavior of P{At)±e 
can be seen in Fig. [51 

4.1. Analysis 

In accordan ce with upper li mit s howed in Fig. [T] and 
results from Ide Souza et al.l (j201l[ ). we expect between 
^ 10^ — 5 X lO'^ events per year in all the sky. The un- 
certainties come from our poor understanding about the 
efficiency with which gas is converted into stars and GRBs 
arc triggered (two unknown factors for Pop III stars). For 
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Fig. 8. Probability for a transient event with duration At 
to be observed by Gaia. At is the time the event stays 
brighter than the Gaia limiting magnitude during the 5 
years nominal mission. 

a good statistics, we create a mock sample of 10^ events 
randomly generated by the Monte Carlo method in order 
to infer the PDF of an event to stay below G = 20 over 
Af(days). The average behavior for on-axis and off-axis 
afterglows as a function of Eiso distribution is shown in 
Figs. [MTOl Once we have P{At), we can generate a sam- 
ple with 10''' events several times and test against their 
probability of being observed by Gaia, given by Fig. [S] 
Combining Figs. |51 [SI and[TUl we obtain the following re- 
sults for the average number of events observed during the 
five years of the Gaia mission: 

- E,,o s; 5 X 10^" 
on-axis:1.34 ± 0.62, 
off-axis: 1.26 ± 0.53, 

on-axis: 2.78 ± 1.41, 
off-axis: 2.28 ± 0.88 

Despite the fact that the total number of on-axis is al- 
ways much lower than the number of off-axis, the observed 
number depends on assumptions regarding the GRB lu- 
minosity functions. For lower energies, the decrease in flux 
due to the observation angle leads to a larger number of 
light curves below the observational threshold. Thus, those 
on-axis have higher probability to be detected than those 
off- axis. 



Eiso < 5x10- 




At(days) 

Fig. 9. PDF of At{days). Probability of an orphan after- 
glow with Eiso ^ 5 X 10^^ ergs to appear above the Gaia 
flux limit for a given time interval. 




At(days) 



5. Discussion 

Despite recent developments in theoretical studies on the 
formation of early generation of stars, there are no direct 
observations of Pop III stars yet. Following the suggestion 
that massive Pop III stars could trigger collapsar GRBs, 
we investigated the possibility to observe their OAs. We 
used previous results from the literature to estimate the 
SFR for Pop III. 2 stars, including all relevant feedback 
effects: photo-dissociation, rcionization, and metal enrich- 
ment. 



Fig. 10. PDF of At{days). Probability of an orphan af- 
terglow with Eiso ^ 10^^ ergs to appear above the Gaia 
flux limit for a given time interval. 

Since we expect a larger number of OAs than on-axis 
GRBs, we estimated the possibility to observe such events 
during the five nominal operational years of the Gaia mis- 
sion. The average number of events observed can be as 
high as to 2.28 ± 0.88 off-axis afterglows and 2.78 ± 1.41 
for on-axis ones. This implies that among the possible 
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afterglows observed by Gaia (|japeli &: Gomboc 2011 ). a 
nonnegligible percentage (~ 10%) might belong to Pop III 
stars. 

However, the detection of those events among the Gaia 
data will not be easy. Gaia will observe more than one 
billion objects all over the sky, and each object will be 
independently detected around eighty times during the 
mission, comprising a total of around 10^^ astrometric, 
spectrophotometric, and spectroscopic observations (after 
detection, the observations are multiplexed in the focal 
plane). Gonsequently, finding the OAs events among all 
that data can be a quite challenging task. 

In Gaia data processing, a system called AlertPipe 
is being implemented to deal with alerts of transient 
events. It is foreseen to operate as follows: first candi- 
date alerts are classified using Gaia data, then sources are 
cross-matched with available catalogues through Virtual 
Observatory or local copies and further classified, and fi- 
nally the alerts are stored on a n alert server and released 
to th e community as VO Events (|Hodgkin fc Wvrzvkowski 
I2OIII : Ivan Leeuwen et ah 201lh . Algorithms are under 
analysis for dealing with GRBs (Wyrzykowski, 2012 priv. 
comm.), but no performance figures are available at the 
present moment. 

Based on Gaia data, the duration of the OA can be 
roughly estimated from the flux variation between two 
subsequent observations: if the event is detected during 
the transit of the flrst telescope, it will be re-observed 
106.5 minutes later when Gala's second telescope re- 
observes the field. Moreover, the light curve will be sam- 
pled several times during the transit of each telescope, 
since at each column of Gala's focal plane an indepen- 
dent magnitude measurement will be performed (measure- 
ments are spaced by 4.4s). 

The light curve alone may not be enough to distin- 
guish between GRB afterglow s and othe r optic al transient 
sources, as noted bv .Japeli fc Gombod (|201l[ ). However, 
as these events have power law like SEDs and no qui- 
escent counterpart, this information should also be con- 
sidered. Further analysis of Gala's BP/RP low-dispersion 
spectrophotomctrjH are needed to distinguish between dif- 
ferent transients with similar characteristics. To perform 
transient event classification, AlertPipe uses several al- 
gorithms, including bayesian cl assifiers, template match- 
ing, and self-organizing maps (jHodgkin fc Wvrzvkowski 
2010h . 



A possible way to search for such objects within a 
large survey is to look for signatures of afterglows from 
Pop HI stars. Two important characteristics of these ob- 
jects are the total energy of Pop HI GRBs, which can be 
much higher than those of Pop I/II GRBs, and the ac- 
tive duration time of their jet, which can be much longer 
than Pop I/II GRB jets due to the larger progenitor star. 
Consequently, the detection of GRBs with very high Eiso 



and v ery long duration c ould be ind i cative of such ob- 
jects (|Suwa fc Iokall201ll: llbma et al.ll201l[) . Thus, they 
should appear as quasi-steady point sources in the radio 
survey observations. But the indication should be com- 
plemented with the constraint on the metal abundances 
in the surrounding medium with high-resolution IR and 
X-ray spectroscopy. Since we do not have any observa- 
tion of these objects, we have to rely on theoretical mod- 
els to compare the data. A way to look for such ob- 
jects that is worth future investigation is the use of auto- 
matic light curve classifiers, which are widely implemented 
for classifying supernoy a e and transients in general 
(IJohnson fc Crottd I2OO6I: iKuznetsova fc Connolivl 12007 



Poznanski et al.ll2007l: iRodnev fc Tonrvll2009t iFalck et al 



2010 : lNewling et al.ll2011:lRichards et al.ll2011l:lSako et al 



2OIII : llshida fc de Souzall2012l) . In principle, the theoret- 
ical model could work as a training set for the classifier, 
which would be then applied to surveys to identify possi- 
ble candidates for further spectroscopical follow up. 

These OAs event will be detected by the Gaia data 
processing pipeline just like any other transient. The 
timescales for raising the alerts are very dependent on 
specificities of the Gaia dataflow, but it is foreseen that, 
in the worst case, the data will be available for analy- 
sis by AlertPipe 24h after the observation. The alerts will 
thus be raised no later than 48h after the Gaia observation 
( Wvrzvkowski fc Hodgkinll2011 ). Nonetheless, it is not yet 
clear if AlertPipe by itself will be able to determine the 
nature of the transient as an OA. 

Moreover, due to the design of the mission dataflow, 
real-time identiflcation will not be possible, but further 
identification of OAs using data from satellites/telescopes 
operating on other wavelengths may be possible, as 
VOEvents will be created by the Gaia data processing 
alert system. Also, OAs could be identifi ed if they trigger 
X-ray detectors, suc h as Swift's BAT ([Barthelmv et al 



20051 ). Fermi's LAT (jAtwood et al.ll2009D . which is fore- 
seen to op erate until 2018, or future instruments, such 



as SVOM ( Schanne et al. 2010h . Finally, the same may 



* BP/RP are Gaia spectrophotometers. BP works between 
330-680nm with 4-32 nm/pixel and RP works between 640- 
lOOOnm with 7-15 nm/pixel. 



also be observed b y other large-scale optical surveys on 
Earth, e.g., LSST '(ilyezic_elalJ l2Q^ and Pan-STARRS 
( Kaiser et al.ll2002 ). improving the sampling of the events' 
light curve and providing information on other optical 
bands. 



6. Conclusion 

It is important to emphasize that our knowledge concern- 
ing first stars and their GRBs is still quite incomplete. 
Many of their properties (e.g., characteristic mass, SFR 
and efficiency to trigger GRBs) are still very uncertain, 
and more reliable informatio n can only come o nce a detec- 
tion is confirmed. Recently, iHosokawa et al. (2011), per- 
forming state-of-the-art radiation-hydrodynamics simula- 
tions, showed that the typical mass of primordial stars 
could be ~ 43M0, i.e., less massive than originally ex- 
pected by theoretical models. Their results, though, are 
affected by assumptions on the initial conditions. This 
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confirms that we are far away from understanding all char- 
acteristics of these objects and any observation would be 
of paramount importance to improve theoretical models. 

In this work, wc estimated the average number of OAs 
events originating from Pop III stars that the Gaia mis- 
sion may observe to be up to 2.28 ± 0.88 off-axis after- 
glows and 2.78 ± 1.41 on-axis ones. In case such events 
are found among Gaia data, valuable physical properties 
associated with the primordial stars of our Universe and 
their environment could be constrained. 
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